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1. Introduction 
Consider the Laplace operator A on the disk D, 

A = dl+d 2 X2 ; D={xeM. 2 :x 2 + x 2 <l} 
and let Uj be the (real-valued, normalized) Dirichlet eigenfunctions 

Auj = —XjUj,x G D; Uj(x) = 0, a; G dD; / Uj(x) 2 dx = 1 

Jd 

with eigenvalues < Ai < A2 < A3 < • • • , indexed with multiplicity. The wave 
trace is the sum 

00 

3=1 

which converges in the sense of distributions. The purpose of this article is to 
announce the following theorem. 

Theorem 1.1. h{t), an infinitely differ entiable junction on 2ix < t < 8, has a finite 
limit and is infinitely differ entiable at t — 2n from the right. 

The significance of 2tt is that it is a cluster point of the length spectrum from 
the left (t < 2tt), as described in more detail below. It is easy to verify that there 
are no geodesies of length in between 2tt and 8; it follows then from |GM] that h(t) 
is smooth in 2tt < t < 8. The content of the theorem is that h is smooth from the 
right up to the endpoint 2ir. The same proof applies to every cluster point 2n£ of 
the length spectrum of geodesic flow on the disk. 

We recall now the relationship between h(t) and the wave equation. Consider 
the initial value problem for the wave equation, 

(d 2 - A)u{t, x)=0, x G D, t > 
u(t,x) =0, x G dD, t>0 
u(0,x) = f(x), x G D 
d t u(0, x) — g(x), x G D 



The solution is 



u(t,x)= / K[ 1 \x 1 y)f{y)dy+ / K<?\x, y)g(y)dy 

J D J D 
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where 



OO 



E 



oo 



Uj{x)uj{y) 




and the trace of the operator with kernel K t is 




oo 



Kt(x, x)dx = ^ cos(^/\jt) 



Re h(t) 



3=1 



The same proof that shows that h(t) is smooth as t —> (2n) + also shows that the 
trace 



is smooth as t — > (2tt) + . 

The close connection between between the length spectrum and the singularities 
of hit) was discovered by way of spectral and inverse spectral problems. In |CdV] . 
Colin de Verdiere showed that on a generic compact Riemannian manifold without 
boundary, the length spectrum is determined by the spectrum (list of eigenvalues of 
the Laplacian with multiplicity) . Duistermaat and Guillcmin |DG] and Chazarain 
[C] showed that the singular support of h(t) is contained in the length spectrum 
and that the two sets are equal generically. Indeed, in the generic case, the sin- 
gularity of h(t) for t near L, the length of a geodesic, resembles a negative power 
of t — L, or, more precisely, a conormal distribution. In case the boundary is non- 
empty, Andersson and Melrose |AM) introduced the notion of generalized geodesic 
length spectrum and proved the analogous inclusion for the singular support of h. 
Subsequent work on inverse spectral problems for domains in the plane and, more 
generally, for manifolds with boundary can be found in |GM| ICdV| El IHeZ) . 

The periodic geodesies on the disk (with reflecting boundary) have lengths 



with k the number of segments (or reflections) of the trajectory and t the winding 
number of the trajectory around the origin. The degenerate cases, (2,1), (4,2), 
(6, 3), . . . , correspond to the trajectory that traverses a diameter 2, 4, 6, . . . times. 
It follows from |GM] that hit) is singular at t = ±Lk : i and smooth in the com- 
plement of the closure of these points. In particular, h(t) is singular when t is the 
circumference of each regular polygon, 



and h(t) is smooth in 2tt < t < 8, since the shortest periodic geodesic with length 
greater than 27r has length 8 = L4 2 (the 2-gon traced twice). 

Thus the issue addressed here that is not addressed in previous works, is the 
behavior of h(t) near a cluster point of the length spectrum. Although we examine 
only the case of the disk, which is far from generic, we expect the analogue of 
Theorem 11.11 to be valid for any convex domain in place of the disk. See the final 
remarks, below. 

It is reasonable to conjecture that hit) has some power law (classical conormal) 
behavior as t — > (2tt) + . What is surprising is that the power is and hit) is in- 
finitely differentiable. The spikes at t = L^.i are asymmetrical and decay rapidly 
for t > Lk,x, so rapidly that an infinite sum of them with singularities at points 




L k/ = fc(2sin(7r£/fc), £ = ±1, ±2,..., k = 2, 3, 4, . . . 



£2,1 < £3,1 < £4,1 < • • • < Lk,i — > 27r k — > 00 
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closer and closer to 2tt still converges at t = 2n . This is explained on a technical 
level by the fact that h{t) is represented by sum of oscillatory integrals in which 
the first derivative of the phase function tends to infinity along one ray. Because 
the derivative of the phase is large, the phase changes quickly, and the correspond- 
ing integrals have more cancellation than one would obtain from a classical phase 
function. 

2. Outline of the Proof 

Our theorem is proved from systematic, optimal symbol properties of the zeros 
of Bessel functions or, equivalently, the eigenvalues of the Laplace operator for the 
Dirichlet problem on the disk. 

Let p{m, n) denote the mth zero of the nth Bessel function J„, p(l, n) < p(2, n) < 
■ ■ ■ . In polar coordinates x\ = rcos6?, Xi = rsin#, the eigenfunctions on the disk 
have the form 

J n (p(m, n)r) cos n9, J n {p{m, n)r) sin n6 

with eigenvalue —A = —p(m,n) 2 . Since the multiplicity of the eigenvalue is two 
when n > 1 and one when n = 0, 



h(t) = V ViM^2(n)e lt "(" l <™> 



E 



for any smooth cut-off functions ipi and 1^2 satisfying 

( , f0, n < -1/4 

, . . 1 1, m> 7/8 , . . , 1 

Mm) = < ' '. ; Mn) = < 1, -1/8 < n < 1/8 

0, m < 3/4 
< ' ' [2, 3/4 <n 

Below we will extend p(m, n) in a natural way to be defined for real numbers m. 
(The extension to real numbers n will be the standard one for Bessel functions.) 
The Poisson summation formula then yields 

h(t) = h kAt) 

where 



Our main result, Theorem 1 1.1[ follows immediately from 
Theorem 2.1. 

<C NltN2 (l + \k\ + \£\)- N2 



d Nl 



dt N ^ 

for 2ir < t < 2?r + 1/10 

The domain of integration for hk,t is the quadrant m > 3/4, n > — 1/4. We will 
deduce the estimates on hk,e from symbol estimates for p(m,n). 

In the sector range of the parameters, m > cn, for any fixed c > 0, the zeros of 
the Bessel functions satisfy ordinary symbol estimates as follows. 

Proposition 2.2. Fix a constant Co > 0. If m > 3/4, n > —1/4 and m > con, 
then 

a) |aj,a*p(m,n)|<C J , fc (m + n) 1 -'-* 
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b) (dmPj d n p) = (tt/ sin a, ctj sin a) + 0((m + n) ) as (m, n) — > oo where a is 
defined by^ 

tan a — a — irm/n 

In the range of values of (m, n) complementary to Proposition ^. 2[ the appropri- 
ate symbol-type estimates involve fractional powers of m and n, and d m p does tend 
to infinity. (For m fixed n — > oo it turns out that <9 m p ~ n 1 / 3 . For our purposes 
the subtlest and most important bound will be the lower bound on the size of d n p.) 

Theorem 2.3. There is an absolute constant c\ > such that if 3/4 < m < cqu 
(co from Proposition \2.2\) then 

a) \diX(p(m,n)-n)\ < C^ k m 2 '^{m + n) 1 '^ 

b) d n {p{m,n) -n)> cim 2/3 n~ 2/3 

Proposition 12.21 and Theorem 12.31 will be proved using asymptotic expansions of 
Bessel functions, which we derive by the method of Watson (jW], p. 251) starting 
from the Debye contour integral representation. 

1 /•00+7T? 

(2.1) JJx)+iYJx) = — / e xsinhz - uz dz 

m J-oc 

We cannot merely quote Watson's asymptotic expansion because we need to differ- 
entiate it. To some extent these differentiated estimates were carried out already 
by Ionescu and Jerison [LJ], but we need quite a bit more detailed asymptotics, es- 
pecially in the transition region where x — nm Cn 1 ^. Also, one needs to choose the 
right coordinate system since differentiation in some directions behaves differently 
from others. 

Theorem 12.11 is proved by integration by parts. Consider the phase function of 
hk,l{t), 

Q = i(tp(m, n) — lixkm — 2ir£n) 

Since Q is smooth, the only issue is the asymptotic behavior as (m, n) tends to 
infinity. It is easy to show from Proposition ^. 21 (b) that for large (m, n) the critical 
points of the phase Q occur near t — Lk t i {a — ni/k). 

The rest of the paper is organized as follows. In Section 3, we carry out the 
proof of Theorem 12. 1[ dividing the (m, n) quadrant of integration into two sectors. 
In the sector where m — > oo, Proposition l2.2l and standard integration by parts and 
non-stationary phase methods apply. In the sector where n — > co, 3/4 < m < CqU, 
the lower bound on d n Q given by Theorem 12.31 is used when integrating by parts 
with respect to n. On the other hand, when integrating by parts in the m variable, 
we will use the oscillation of e tkm only and include the rest of the factors in the 
exponential in the amplitude. 

In Section 4, we prove the symbol estimates and asymptotic expansions for Bessel 
function following the method of steepest descent (Debye contours). In Section 5, 
we deduce the symbol estimates for the zeros of Bessel functions stated above. 
We conclude with remarks about the relationship with earlier work and about the 
methods that can be expected to lead to the analogous result when the disk is 
replaced by a convex domain. 



We extend the definition of a continuously across n = 0, a = n/2 by the reciprocal equation 
l/(tana — a) = n/nm. 

4 



3. Integration by Parts 

We now deduce our main estimate, Theorem 12.11 from the symbol estimates for 
p(m,n), Proposition 12.21 and Theorem 12.31 Denote 

Q = i(tp(m, n) — 2ixkm — 2ittri) 

Then 

(3.1) e« = - =1- d n e* 

to n p — 2n£ 

(3.2) e Q = e i(*P-2TtoJ_ 9me -27rim& 

2nk 

We will divide the region of integration into sectors. Consider first the "non- 
classical" region 3/4 < m < cqu. Define 

Z w (f)= f ^ 1 {m)^ 1 {cn/m)e i ^ p{m ' n ^- 2wkm - 2 ^dmdn 
Jr 2 

We focus at first on the case I = 1, Q = i(tp(m, n) — 2irkm — 2irn) from which all 
the singularities near t = 2ir arise. Applying formula [3. II and integrating by parts, 

4,1 (t) = / ipiim^x^n/m)— % ——d n e Q dmdn 



td n p — 2n 

= f i — r ^^^i{m)^2{cn/m)e'^dmdn 

J {tdnp-2-KY 



+ / — t/^! (m)-0^ (cn / m) — —e^dmdn 

m to n p — 2tt 



Repeating we obtain 



ifc,i(*) =c l / ^^i{m)^"(cn/m)-— —-re Q dmdn 

' ~~ z [to„p — 2-k) z 

C2 / — ipi(m)ip[(cn/m)— - — ""^ e^dmdn 
m (to n p — 2 r K) 1 



+ / %l)i{m)ipi(cn/m) 



(td nP - 2-kY (td nP - 2tt) 3 



e^dm dn 



for appropriate coefficients Cj(t) (polynomial in t). After N integrations by parts, 
Ik,i is expressed as a linear combination of terms with integrand 

6 (^-2^^.-^ i (« 2 + 2«3 + • • • + Na N+1 = N) 

times cutoff functions ipi(m)i/)i(cn/m) or derivatives of these cutoff functions. For 
t > 27r, the denominator has the lower bound 

td n p -2tt> 2n(d n p - 1) > cm 2 / 3 n" 2 / 3 

from Theorem 12. 3f b). Moreover, Theorem 12. 3f a) says, in particular, that 

Kp\ < m^n 1 / 3 -! 
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Denote B = ai + 03 -| — ■ + Ozv+i- Then each of the integrands is bounded by 

(^ri a2 Kp)" 3 '"(3: + vr +i 



(td n p - 27r) Ar + a 2+ a 3+---+ Q «+i 



(j7j3 n~ 3 ^A r +a 2 +a 3 H ha N + l 

m f B n iB n -(2a 2 +3a 3 + --- + (A'+l)a J v + 1 ) 

^ m 2/3 n -2/3^7V TO |-B n - JS 
n -(a 2 +2a 3 +---+A r a N+ i) ^-JV 



( m 2/3 n -2/3)iV ( TO 2/3„-2/3)iV 

— 2/V -iff 

= to 3 71 3 

In particular, Ik,i(t) is represented by a convergent integral. 

Next, to prove rapid decay in fc we integrate by parts in to using substitution 
O The first step is 

h,i{t) = J e^ 2 ^^d m e- 2 ^ mk Mm)Mcn/m)dmdn 

nl P e Q-ip 1 ( m \ip 1 (cn/'m)d'm dn 
2nk 

e^^\(rn)^>i(cn/m)drndn + / —e^^i(m)^\(cn/m)dmdn 

2nk J linn k 

The term with the factor n/m 2 also has ip[(cn/m) so that it is supported where 
ctl ~ to and the term n/m 2 is comparable to 1/n ~ 1/m. All terms have a gain of 
a factor 1/fc except the one in which the derivative d m falls on p. In that case, 

[SW.Pl 3 ^to' 1 / 3 
In all, one step of type (|3.2I) yields the factor 

n l/3 m -l/3 

k 

Now we consider systematically what happens when steps of type (|3.2I) are applied 
after N steps of type (|3.ip . If the derivative J m in the integration by parts falls on 
a factor d J n p, then this gets replaced by d J n d m p/k and thus the bound is improved 
by the very favorable factor 

1 

mk 

If the derivative d m falls on the e l ('p~ 27rn ) as i n the first step, we have, as before a 
factor 



k 

If the derivative falls on a cutoff, then the gain is 1/fc. Finally, if the derivative falls 
on the denominator (td n r — 2tt) then it produces a factor 

1 d n d m p , m _1 / 3 n" 2 / 3 _ 1 



k{td n p — 2n) ~ fcm 2 / 3 n~ 2 / 3 fcm 

In all, the worst case is the factor n 1 ^ 3 m~ 1 / 3 /fc for each integration by parts in to. 
Thus if we integrate by parts N times in n and M times in to, the integrand will 
be bounded by 

m-i N n-i N ni M m-i M k- M 
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Therefore, if we choose N sufficiently large that 

\N-\M>2 

then we obtain a convergent integrand that gives a bound on the integral by k~ M . 

For I ^ 1, the bound is much simpler. d n p is very close to 1 for small m/n, so 
if 2tt < t < 4tt — S for any fixed S > 0, then the denominator in the integrands, 

\td n p - 2n£\ » \£\ 

for all integers i ^ 1. By integrating by parts AT times, one finds the bound 

for each N. The rapid decrease in k follows from similar reasoning to that given 
above for 1^,1- 

Next, consider derivatives {d/ dt) Nl Ik.i- The integral representing this expres- 
sion just has an extra factor of p Nl in the integrand. This extra factor has size 
(n 1 / 3 ™ 2 / 3 )^ and symbol type bounds of the obvious kind after differentiation with 
respect to m and n. Thus one can compensate for these higher powers by more 
integrations by parts, and nearly the same proof as above shows that the derivatives 
of Ik,i(t) are also rapidly decreasing in (k,£). This ends the main portion of the 
proof. 

What remains is to make estimates for the integrand in the region region m > cn 
of integration. 

(3.3) hk.e — 2Ik,e = /h/'i( m )V'2(?'0 — 2ty\(m)^\(cn/m)\e®dmdn 

For this region we use Proposition 12.21 The informal idea is as follows. Fix a > 0, 
and consider a ray in (to, n) space defined by 

tan a — a ~ nm/n 

If VQ — > (0, 0) as (to, n) — ¥ oo along this ray, then the asymptotic formula implies 

tn/ sin a — 2Trk; iasina = 2ir£, 

Thus if Q has a "critical point near infinity" we can solve these equations for a and 
t and find 

a = n£/k; t = 2k &in(2ir£/k) = L M 

This explains the singularities at t = L^^. The fact that the phase is nonstationary 
at other values of t will lead to a proof that hk.i{t) is smooth at each point t ^ L^j. 

In more detail, first consider the range |(m.,n)| < C, truncating the integrand 
with a smooth bump function in (to, n) variables. In that range, since the deriva- 
tives of p are bounded, 

\d n Q\>c\l\, \d m Q\>c\kl 

for sufficiently large \k\ and \£\. Hence, writing 

e Q = (l/d n Q)d n e Q , e Q = {l/d m Q)d m e®, 

and integrating by parts, one finds that the integral decays like l/(|fc| + \£\). Re- 
peating N times, one finds that the integral decays like 0(l/(|fc| + 1^1)^) for any 
N. 
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Next we turn to the range |(m, n)| > C, m > cn. In this case, we will also do 
integration by parts either in the variable m or n, and use lower bounds on |<9 n <2| 
or \d m Q\. 

Here we will use the asymptotic formula for (d rn p, d n p) of Proposition 12 . 21 b and 
restrict t to 2ir < t < 2ir + 5 for suitable small number 5. 
First we confirm 

(3.4) \d„Q\ = |ta/sina- 2n£ + 0(l/(m + n))\ > c(\£\ + 1), 

Recall that a is defined by tana — a = irm/n. If n > 0, then, since m > cn, we 
have (Xq < a < t:/2 for some fixed ao > depending on a In the remaining range, 
> n > — 1/4 and m > C — 1/4, which implies that it/ 2 < a < n/2 + S for some 
small 6 of size on the order of 1/C. It follows that 

1 + ao/10 < a/ sina < vr/2 + 45 

The lower bouncQ on a/ sina bounds |9 n Q| from below when £ < 1 and the upper 
bound on a/ sina bounds |<9„Q| from below in the case £ > 2. Thus we have proved 
(|3.4[) . We use this bound and integration by parts in n in the range. \£\ + Ci > \k\. 
Lastly, if k > \£\ + C2, then (taking C% = 10/ sinao) we have 

(3.5) \d m Q\ = |*7r/ sina -2ttA; + 0(1/ (m + n))\ > c(\k\ + 1). 

On the other hand, if k < 0, then (|3.5[) is obvious since <7r/ sinao > 2ir 2 . In the 
range |fc| > \£\ + C2, we use (|3.5[) and integration by parts in the m variable. This 
concludes the proof of Theorem 12.11 

4. ASYMPTOTICS OF BESSEL FUNCTIONS 

In this section, we establish the optimal symbol properties of 

H„(x) := J n (x) + iY„(x) 

as a functions of two variables (x, n). (The function H u (x), known as a Hankel 
function or Bessel function of the third kind, is denoted (x) in Watson's treatise 
|Wj p. 73.) The asymptotic formula for the Bessel functions as the order and 
variable tend to infinity was discovered by Nicholson in 1910. In 1918, Watson used 
the Debye contour representation to give an appropriate bound on the error term. 
In his treatise on Bessel functions ([W] p. 249), Watson says of his own method 
that it is "theoretically simple (though actually it is very laborious)." To prove 
Theorem 12.31 we will carry out the even more laborious process of differentiating 
Watson's asymptotic formulas. 

To state the symbol properties of H n (x) in the sector n < x < 2n, especially in 
the so-called transition region in which x is very close to n, will require a different 
coordinate system (/3, v). For x > n > 0, we write v = n and define j3 by x cos j3 = n. 
Define a(v, (3) by 

(4.1) R v {y sec/3) = e lt/(tan ' 3 ^' 3) a(^, 0) 

Proposition 4.1. If v > 1/2 and < /3 < n/4, and a(p,/3) is defined by (|4.ip . 
then 

2 The key here is that in the case 1=1, tas'ma — 2tt > 2^0^/10 > 0. We have avoided the 
critical points associated with t = Lj. 1 — > 2ir. They occur at infinity along a rays in (m, n) space 
near the n axis, rays that are not in the sector m > cn. 
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We will also need more detailed asymptotics involving Airy functions. Define 
the function A(y) for y £ M as the solution to the equation 

A"{y) + 2yA(y) = 

with initial conditions 

A(fi) = r(l/3)6- 2/3 (3 + tV3); A'(0) = r(2/3)6- 1/3 (-3 + iVi) 

Proposition 4.2. Denote b(v,f3) = e - lu{tan3 p)/3 v^ 1/3 A(y) withy = (l/2)v 2/3 tan 2 (3 
and a(y, /3) from Proposition \4-l\ Then for < f3 < 7r/4, 

a) \a- b\ ^ i/- 1 . 

b) \d v {a-b)\ <v- 2 . 

c) \d p [a-b)\ ^ v-' 2/3 [also < l/vf3 if f3 > v- 1 ' 3 ] 

Corollary 4.3. If v » 1 and < (3 < n/4, and a(u, f3) is defined by (|4.ip , then 

a) \a\ w v~ 1/2 f3~ 1/2 provided i/" 1 / 3 < (3 < it/ A. 

b) \a\ w i/" 1 / 3 provided j3 < v~ 1 / 3 . 

In the remainder of the section, we will prove Propositions 14.11 and 14.21 and the 
corollary. The range, x > 2n will be discussed at the end of the section. 
Proof of Proposition^. 1\ Define the phase function ip(z) by 



tp(z) = — (irsinhz — vz) = sec/3sinhz — z 
v 

The contour integral (|2.ip can then be written 



-1 fOC-fK'l 

H u (vsecf3) = — e vip(z) dz 



7" J-oo 

The contour of steepest descend passes through z — i/3 and is parametrized by the 
two curves, z = Ci(V, f3) + i/3 and z — (2(1", f3) + i/3, in which the functions Q(r, f3), 
j = 1, 2, solve 

(4.2) (p(( + i[3) - tp(if3) = (sec/3) sinh(C + if3) - C~ itan/3 = -r, r > 0, 

and satisfy Ci(0,/3) = C2(0,/3) = 0, ReCi(r,/3) < 0, ReC2(r,/3) > 0. Moreover, 
(i(r, (3) — > —00 — i/3 and ^(r, (3) — > 00 + i(7r — /3) as r — > 00. Thus 

p Utp(ifi) rca+TTi 

H v {vsecp) = — / ^M')-<pm)dg 



TTI 

e 



—00 

ii/(tan/3-/3) /-oo 



e-^(9 r C2(r,/?)-a r Ci(r,/3))dr, 







717 

and we have derived the formula for a(v, f3), 

(4.3) a(v,[3) = -. I™ e-» r (d r ( 2 (r,f3)-d r C 1 (r,{3))dr 

™ Jo 

Lemma 4.4. For j = 1, 2 , 

' r -i/a j8 -i/2-fc j r < ^3 



r -(2+fc)/3 ) /3 3 < r 



3 We follow IW| p. 244 and pp. 249—252, except that where Watson uses e xr , we use e 
our expressions differ from his by factors x/v = cos/3. 
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Lemma 14.41 follows in a straightforward way from implicit differentiation and 
induction, but the proof takes a few pages. Abbreviate by £ the functions Q(r,f3) 
along the contour of steepest descent satisfying (I4.2[) . Differentiating (|4.2p with 
respect to (3 yields 

[cosh(i/3 + C) - cosh(i/3)]<9,3C + (tan f3) sinh(C + i/3) +i cosh(C + i/3) - i sec (3 = 

Therefore, using cosh(a + b) = cosh a cosh 6 + sinhasinho with a — ( + i(3 and 

b = -i/3, 

— i(coshC — 1) sec/3 
^ = cosh(C + i/3) - cosh(^) 

Define 

Fi(C) = l/sinh(C/2); F 2 (C,/3) = 1/ sinh(C/2+i/3); G(C, /3) = (coshC-l)(sec/3)/C 2 
Since 

cosh(^ + C) - cosh(i/3) = 2 sinh(C/2) sinh(i/3 + C/2), 

we may write 

(4.4) fyC - y FiF2G(2 
Similarly, 

(4.5) d r ^ = ^-F x F 2 cosp 

Lemma 4.5. Along the contour of steepest descent C = Ci(r, /3) and £ = £2(^1 P), 
for all k>0, 

'ICI-'GS + ICI)-* for |CI<1 
e-l^l/ 2 /or (CI > 1 

VlSl»«r M <{\W + ^- k -' lor KjSl 
[1 /or ICI > 1 

To begin the proof of Lemma |4~5| denote by D^F{C, : /3) the derivative of F((, /3) 
with (3 fixed and let DpF((,f3) represent the derivative of F(£,/3) with C fixed. 
(This is to distinguish from the partial derivative dpF(((r, (3), (3) representing the 
derivative with r fixed.) We bound and Dp derivatives of F\, F2 and G as 
follows. On the curve of steepest descent (C = Q(r, (3), j = 1, 2) 



-< 



(4.6) |£>pi(C)| = |(d/dC)^i(C)l 



-< 



icr fc_1 for ici < 1 

e -KI/2 for |^| > ! 



(4.7) |^F 2 (C,^)U| e _ |c|/2 for |C| >l 



(4.8) |D^G(C,/3)|^ 



1 for |C| < 1 

^eKI/ICI 2 for |C| > 1 
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To prove (|4.6[) when p = 0, note that if £ = £ + 297 and £ > 0, then 

| S inh( C /2)|>(^-l)/2>|^ /4; ^° 

Along the contour C = £ + if) = ( 2 {r, /3), < 77 < tt - /3. (Recall that < /3 < tt/2.) 
Furthermore, the slope of i] as a function of £ is at most V3 (jW] 8.32, p. 240) so 
that C > ICI - tt and £ > ICI/10- It follows that \F X (C)\ = |l/sinh(C)| < ^ICI" 1 for 
all C = (2(r,/3). For |C| > 20, we have £ > 2 and consequently \Fi(z)\ < 4e~ 5/2 < 
4e ^/2 e -|CI/2^ Thig proveg gjj) for p = o and ^ = The Qther branch £ = ^ 

is similar, with the only difference that the curve Ci is in a horizontal strip of the 
complex plane is below the £-axis: — j3 < 77 = Im £i < 0. The case of p = 1,2,... 
are easy consequences of the same estimates. A very similar proof to the one for 
PTSjl gives (gTTj) because 

|i/3 + C|«/3 + |C| 
along the contour. The estimate (I4.8[) is easy. 

We can now carry out the proof of Lemma 14.51 by induction. Parts (a)-(c) for 
k = are the same as pjjjflTgf for p = 0. Part (d) for k = (k + 1 = 1) follows 
from (a)-(c) for fc = and the formula (I4.4p . 

For the induction step, assume (a)-(d) are valid up to k. Part (a), <9o +1 .Fi(£(-r, /?)), 
is a linear combination of terms 

(D^dfOidfO ■ ■ ■ {dfO (91 + • • • + q P = k + 1; ft > 1) 

Part (b), 9^ +1 J 7 2(C( r j ft), ft), is a linear combination of terms 

(D P c D s fi F 2 (z, P)){dfQ{dfQ ■ ■ • (^ P C) (ft. + • • • + ft + a = fc + 1; ft > 1) 
Part (c), d^ +1 G(C,{r, (3), (3), is a linear combination of terms 

(D P c D s ? G(z,0))(dfQ(dfO ■ ■ ■ (8fQ ( qi + ■■■ + q p + s = k + 1; y > 1) 

Since the induction hypothesis says that d%C, has the given bounds for all 1 < q < 
k + 1, each of the factors has appropriate bounds and multiplying them yields the 
appropriate bounds for (a)— (c) with k replaced by k + 1. 
Lastly, to prove the induction step for (d), observe that 

dp l (dpQ = dp l [(-i/2)F 1 F 3 GC 2 ] 

is a linear combination of terms of the form 

{df Ft (0) {df F 2 (CIS)) [df G) {df C) [df C) 

with ki + k 2 + • • ■ + fcs = k + 1. Since fcj < fc + 1, we have already proved the 
appropriate bounds on each factor. Multiplying them yields the correct bound for 

ak-\ " 



dg +2 C- This ends the proof of Lemma T4.5I 



To convert the implicit bounds of Lemma [4. 5 1 to ones in terms of r and /3, observe 
that for £ = Cj( r i 0)i 

r r i/2 j9 -i/ 2j r <p 

ICI « < r 1 / 3 , /3 3 < r < 1 

I 1 + log r, 1 < r 

In the range r > 1, there is also a more precise estimate along the contour, namely, 

e ici ~ e i?i ~ r (C = £ + 
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It follows that 

r < /3 3 

/.>' + (, <-k - { /-'■ ; . /3 3 <r < 1 

1 < r 

With these upper and lower bounds on |^| and /3 + |£|, one can rewrite Lemma 
1431 as 




Lemma 4.6. Along the contour of steepest descent £ = £i(r,/3) and £ = C2 ? 
for all k>0, 

' r -l/2ol/2-As 

a; |^§Ji(C(r,jS))| ^ ^ r-iC^ 1 ) 



6; |^F a (C(r, /?),/?) | 
c) |^ fe G(C(r,/3),/3)| 




' r l/2p-3/2-k r < fiZ 

d)\dl +1 C,{r^)\< lr-^ k /3 3 < r < 1 

1 1< r 



It is now routine to confirm Lemma 14.41 Differentiate (14.51) , 
dp r C = (-l/2)d k (F 1 F 2 cosp) 
The right-hand side is a linear combination of terms of the form 
{dfF 1 ){dfF 2 ){df cos/3), k x + k 2 + k 3 = k 



which are bounded using Lemma 14.61 Now that Lemma 14.41 is proved, Proposition 
14.11 follows from the formula (|4.3[) for a(u,/3). 

We now turn to the proof of Proposition ^. 21 The function (p from (14. 2[) can be 
rewritten 

(4.9) (p(z + - (p(if3) = (i tan /}) (cosh z-l) + (sinh 2-2) 

The Taylor approximations coshz -1» z 2 /2 and sinhz — z — z 3 /6 so the cubic 



(4.10) $(z) = {itan/3)z 2 /2 + z 3 /6 

Following Watson, we prove Proposition ^. 21 by comparing £1 (r, 0) and (, 2 {r, /3), the 
two solutions to (|4.2p , to the functions Zi (r, /?) and Z 2 (r, /3) solving the correspond- 
ing equation 

(4.11) $(Z) = -r, r > 

with Zi(0,/3) = Z 2 (0,/3) = 0, and ReZ x (r,/3) < 0, ReZ 2 (r,/3) > 0. 
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Following Watson again, curve Z\(t) — > — oo — itan/3 as r — > oo, whereas Z<i(r) 
is asymptotic to the ray whose argument is 7r/3 as r — > oo. Thus we define 

e v * (z) dZ 

-oo— i tan /3 

In order to show that this function b(/3, v) is the same as the function b of Proposi- 
tion ^^! we evaluate it using the contour consisting of the two rays Z = —i tan/? — £ 
and Z = — itan/3 + £e I7r / 3 , £ > 0, rather than the steepest descent contour defined 
using Z\ and Z-2- With the parametrization given for these rays, and the change of 
variable £ = u/v 1 ^ , one obtains 



7e 7(^tan 2 /3)/2 + £ - (u£ tan 2 /3)/2 



(4.12) 603,1/) = e -^ tan J «/ 3 / e -^ J /6 

./0 

= e- i ( l ' tan3 «/ 3 l /- 1 / 3 ^((^/ 3 tan 2 /3)/2) ( 7 = e"/ 3 ) 



where A(i) is defined by 



A(t) 



We claim that 

Indeed, for A € C, define 



e -« /6[ 7e 7*« + e -ti 
A"(i) + 2tA(t) = 

/>OC 

F x (t) = / e - u3/6+xtu du 
Jo 



]du 



Then 



Hence, 



-1 = / {d/du)e- u / 6+xtu du 

{-u 2 /2 + \t)e- u3/6+xtu du 
= (-l/2\ 2 )F{{t) + \tF x {t) 

--^) + A 3 iF A (0=-A 2 



■1 give the equation A"(t) + 2tA(t) = 0, 



Then A(t) = jF 7 (t) + and 7 3 

as desired. 

Thus A(t) is an Airy-type function, identified uniquely by its value and derivative 
A(0) - r(l/3)6- 2/3 (3 + iV3); A'(0) - r(2/3)6~ 1/3 (-3 + iVS) 

Writing A in terms of its real and imaginary parts, A(t) = u(t) +iv(t), we find that 
the Wronskian takes the constant value 

(4.13) u(t)v'{t) - u'(t)v(t) = u(0)w'(0) - u'(0)u(0) = T(l/3)r(2/3)V3 = 2tt 

Having identified b with the function of Proposition 14.21 we can now proceed 
with the proof. 

Lemma 4.7. For j = l, 2, < f3 < ir/4, 

r i/3 r > fi 3 
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\Zj\< 



Lemma 14.71 is routine and the proof is omitted. In fact Q grows more slowly than 
Zj for for r >> 1 (like logr), but we don't make use of this. 

Lemma 4.8. For j = 1, 2, 

f ,.1/201/2 r</? 3 

a) \(p'(Cj)\ = |itan/3sinhCj +coshCj - 1| £ < r 2 ' 3 f3 3 < r < 1 

[ r 1 < r 

&; |$'(^)| = |itan/3Zj +Z 2 /2| £ 

Part (a) of this lemma is proved in [IJ] 9.15 and 9.16, p. 1072. The proof of part 
(b) is similar and is omitted. 

Lemma 4.9. For j = 1, 2, < j3 < tt/4, 

! r 3/2^-3/2 r < ^3 

r /3 3 < r < 1 

f l/3 ^ < r 

Proof. Fix a small number ro > 0. For Or > ro, the estimate follows immediately 
from Lemma l4~7l For this argument we will use subscripts on the constants, because 
some will depend on others. For \z\ < 1, \<p(z) — $(z)| < Ci(/3|z| 4 + |z| 5 ). Fix 
r, j3 and j = 1 or 2. We consider first the case < r < (3 3 . Denote tUj = 
Q(r, (3) — Zj(r, (3) and Zj — Zj(r, (3). Our goal is to show that 

\wj\ < Cr 3 ' 2 ^ 3 ' 2 

The root in the appropriate half-plane of 

ip{w + Zj) + r = 

is w = Wj. We will show that for a suitable constant C to be chosen later, the curve 
S = {ip(w + Zj) + r : \w\ = Cr 3 / 2 ^ 3 / 2 } 

encloses the origin. Thus the root Wj is inside. We will show that S enclosed the 
origin by showing that ip(w + Zj) + r — &(Zj)w is suitably small. Recall that 
W{Zj)\ > cir 1 / 2 /? 1 / 2 . It follows that the circle 

So = {&(Zj)w : \w\ = Cr 3 ' 2 ?- 3 ' 2 } 

has radius at least Cc\r 2 j 'f3. 

By Lemma S3 \Zj\ < C 2 r 1/2 /T 1/2 . We will require that 

(4.14) Cr 3 / 2 f3- 3/2 < C 2 r l ' 2 r 112 

so that \w + Zj\< 2C 2 r 1/2 p- 1/2 and 

\<p{w + Zj)-${w + Zj)\ < d{/3\w + Z 3 \ 4 + \w + Zj\ 5 ) < C 3 r 2 /f3 

with a constant C3 depending only on C\ and C 2 (using r < j3 3 ). Furthermore, 
since F(Zj) = — r, 

$(w + Zj) =-r + &(Zj)w + (1/2) (i tan £ + Zj)w 2 + w 3 /6 

and 

|(ttan + Zj)w 2 /2 + w 3 /6\ < C 2 f3- 2 r 3 + C 2 C 2 r' 1 ' 2 7/2 + C 3 r 9 / 2 /T 9/2 
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r < (3 3 
r>(3 3 



We will require 

(4.15) C 2 /3~V + C 2 C 2 r 7/2 /r 7/2 + C 3 r 9/2 /T 9/2 < {l/W0) Cl Cr 2 /(3 
Now our final requirement on C is the lower bound 

(4.16) C 3 < (l/4) Cl C 
Combining these estimates, we have for \w\ — Cr 3 / 2 f3~ 3 / 2 , 

\ip.(w + Zj) + r - $'(Zj>| < C 3 r 2 /(3 + (l/100) Cl CV 2 //3 < (l/2) Cl Cr 2 / f3 
so that 

\(p(w + Zj) +r\> (l/2)ciCV 2 //3. 

On the other hand \§'(Zj)w\ > c\Cr 2 / '/3. So 5 is a loop surrounding the origin. 

Finally, we check that all three requirements on C are satisfied. For (|4. 16|) fix 
C = AC 3 /c\. Now that C is fixed, we choose r > sufficiently small that the other 
two inequalities (|4.14|) and (|4.15p are satisfied for all r, < r < r , < r < /3 3 
(0 < /3 < 7r/4). The remaining case, /3 3 < r < tq is similar and slightly simpler. It 
will be omitted. This concludes the proof of Lemma [ 



Lemma 4.10. For j = 1, 2, 

o;i^i<| r _ a/8 r>p3 

b) \d r Q < \r~ 2 / 3 /3 3 <r<l 

[r 1 l<r 

c) \driCj - Zj)\ < 1 oHr 

Proof. Differentiating (|4.1ip . we find 

^(Z^drZi = -1; tfi&drCj = -1 

and the bounds of Lemma 14.81 imply parts (a) and (b) . For part (c) 



$'(^V(0) $'(W(Ci) 

Lemma 14.81 implies 

|$'(^)v'(Ci)|fcmax(rAr 4 / 3 ) 
The formula of $' and Lemmas 14.71 and 14.91 imply 

|$'(Z,) - $'(C,)| ^((3 + \Z S \ + U)\Z 3 — 0| 3 < r 4 / 3 /3 3 < r < 1 

[r 2 / 3 1 < r 

Hence, - 3 1 Lemma gTT] implies 



r 4/3 ^3 < 



Hence, similarly, - V(G)I 3 1- This concludes LemmaEEU] 
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a) \d r df3Zj\ < 

b) \d r dpQ < < 



Lemma 4.11. For j = 1, 2, 

1^-1/2^-3/2 r < ^3 

\ r -l ft 3 < r 

V-l/2^-3/2 r < ^3 

r -i /3 3 < r 

cj |5 r ^(0 - Zj)\ < r- 1 / 3 
Proof. Differentiate the implicit equation with respect to ft to obtain 

dpZj = -(isec 2 ft)Z 2 j /2^'{Z j ) 

and 

dpdrZj = -(isec 2 ft)d r Z j [2Z j ^'(Z j ) - F"(Z j )Z?]/2F'(Z j ) 2 

We have already bounded each of these terms and the bounds combine to give 
Lemma 14.111 (a) . To prove (b) , differentiate (|4.2p with respect to ft to obtain 

dpd = -*(scc 2 ft)(co S hQ - W(Ci) 

Then, differentiating with respect to r, 

d r d b Q = -i(sec 2 ft)d r (j[ S mh( j( p'((j) - ¥>"(<,■) (cosh £ - 1)]M&) 2 

The estimates above for t? r £j-, an d 

!/3 r < /3 3 

r l/3 /33 < r < ! 
r 1 < r 

combine to give part (b) of Lemma 14.111 
To prove (c) write 

(4.17) 

dp(Q - Zj) = -i S ec 2 ft[Z 2 ( V '((j) - V '(Z 3 )) + Z 2 ^'(Z 3 ) - &(Zj))+ 

(4.18) (Zj - Q?&{Zj) - 2 ( cosh G - 1 - C-/2)$'(^)]/*'(^>'(0) 

With z is a point on the line segment from Zj to Q, 

(4.19) 

|0/j(C- - ^)l < C\Zj\ V'WII^ - 01 + 1^| 2 (/3|^-| 3 + 
(4.20) |z,| 4 ) + |z,||^ - + I0l 4 ll^(^-)|]/l$'(^)^(0)l 

Using the preceding bounds and (|4.19[) . 

\»(r 7^/ r3/2 ^ 5/2 r ^P 3 
1^(0-^)1^^/, ^3< r <! 

In the range < r < 1, differentiation of (14. 17)) with respect to r replaces terms Zj 
(with the bound r 1 / 2 /? -1 / 2 ) by (with the bound r~ x l 2 ft~ x l 2 ) or the similar 

replacement of Q with d r Q. This results in a bound of the same type as (|4.19[) 
with an extra factor of 1/r. In other words, 

y/2^-5/2 r < /? 3 

Z,)\ r- {r-V 3 ft 3 <r<l 
- 1 1 < r 
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(The case r > 1 follows separately from parts (a) and (b).) In all three cases this 
is less than r -1 / 3 , so this concludes Lemma [4.111 
Define 



b(p,v) = / e- vr (d r Z 2 {r,p) - d r Z x {r,p))dr 
Jo 

Then Lemmas 14.101 and 14.111 imply that 

poo 

(4.21) \a(fi t v)-b(J3,u)\< e~ vr dr<l/v 



(4.22) \d v [a(fi,v) - b(fi,v)]\ < / re^dr < l/v 2 

Jo 

(4.23) \df>[a(p,u)-b{p,v)]\< r-We-^dr < v~ 2 ' 3 



This concludes the proof of Proposition ^. 21 

Next, in order to deduce Corollary 14.31 we will prove some estimates for b(is,(3) 
that will also be needed in the next section. 



a) \b\ « 

b) \dub\ 3 



Proposition 4.12. If y > and v > 1/2, < /3 < tt/4, then 
'i/-V3 pKv-V 3 

v -l/20-l/2 0> v -l/3 

' V 4 / 3 /3<Z/-V3 
v -l/20S/2 £ > „-l/3 

y/ 2 p 3 / 2 ^>i/-V3 

Proof. The proposition will follow easily from the formula for 6 in terms of v4(y) 
and the estimates for y > 0, 

(4.24) |A(y)|«(l + lO- 1/4 5 K(2/)l«(l+y) 1/4 

(|4.24p is well known, but we include a sketch of a proof. The upper bounds are 
standard; indeed, the asymptotic behavior as y — > oo follows from the fact that A(y) 
is a multiple of the Hankel function H 1 / 3 ((2y) 3 ^ 2 /3) [Wj p. 252. The lower bounds 
(for all y > 0) follow from the upper bounds and the fact that the Wronskian (|4.13[) 
is constant. 

Next, using (|4.24j) we deduce (a). Recall that from (|4.12l) . 

\b(v,p)\ = V-V*\A(V)\ « ^ 1/3 (1 + M)~ 1/4 ; y - (l/2)^ 2/3 tan 2 (S 
When /? < i/- 1 / 3 , \y\ < 1 and |6| w z/- 1 /^. when /3 > i^ 1 / 3 , 

|6(i/,)8)| « ^- 1 / 3 ( I / 2 / 3 /3 2 )" 1 / 4 « v' 1 /*?' 1 /* 

Parts (b) and (c) of Proposition 14.121 are proved as follows. Differentiating (I4.12[) 
gives 

\d v b{v, 0)| 3 /3 3 ^- 1/3 |A(t)| + ^- 2 / 3 2 |A'(t)| < miniis- 1 / 2 ^/ 2 ,^ 3 ) 

and 

3 v 2/3 P 2 \A{t)\ +v 1 /' i p\A'{t)\ < mm^ 1 / 2 /? 3 / 2 ,^ 1 / 3 ^) 
Proposition 14. 121 (a) and (|4.21l) imply Corollary 1431 

Finally, we discuss the range 7r/4 < j3 < tt/2 and beyond. To formulate this we 
return to the variables (x,n). Fix c > and let @o > be the smallest number such 
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that cos/3 = 1/(1 + c )- Let tt/2 < (3i < tt satisfy cos/?i = -1/4. For x > 1, 
n > —1/4 and x > (1 + c)n, define /3(ie, n) as the unique number /?o < /3 < /?i such 
that a: cos f3 = n. 

Proposition 4.13. Let a: > 1, n > — 1/4 and x > (l + c)n for a fixed c > 0. Define 
a(x, n) by 

H n {x) = e l{xsinli - nf}) a{x,n). 
with f3 = (3{x, n) defined above. Then 

\did k n a{x,n)\<x- l l 2 ->- k 

To explain the connection with the previous notation, if v = n and x = v sec f3, 
then a{v,(3) ~ a{x,n). The distinction between the coordinate systems is that 
d v is the derivative with j3 (or equivalently xjv fixed), whereas d n represents the 
derivative with x held fixed. The d v direction is special when f3 is near and the d n 
direction is special when f3 is near 7r/2. The estimates we carried out in the range 
uO < (3 < 7r/4 can be extended to [3 — > tt/2, but they require additional factors of 
sec (3 which tends to infinity. They do not suffice: In the eventual analysis of the 
behavior of zeros p(m, n), the range x = p(m, n) > (1 + c)n corresponds to m > cn 
and estimates in the (v, (3) coordinates give rise to error terms of size 0(l/n) when 
what is needed is 0(l/m) — 0(l/p). 

Proposition 14.131 was already proved in the case n > 0, j = 0, k = 0, 1, 2 in 
Theorem 9.1 (i) of [TT]. (In the notation a x (v) of [D], v — n, a x (n)x^ 1 ^ 4 (x — 
rt) _1//4 = a(x,n), and in the range of variables specified here, x — n « x.) The 
full proof of Proposition 14. 1 31 follows the same procedure as in [LI pp. 1068-1072, 
with the only extra ingredient being the systematic treatment of derivatives of all 
orders, which was already carried out above in the very similar proof of the symbol 
estimates for a(y,j3) in Proposition ^. II (a). These details will be omitted. We call 
attention to one difference. The integrals in the proof of Proposition 14.11 involve 
e~ vr dr as v — > oo, whereas in the proof of Proposition 14.131 (following [LI]) the 
integrals involve e~ xr dr and x — ¥ oo. We mention this in order to explain why the 
proof is unchanged when the range of n is extended from n = i/>0ton> — 1/4. 
The range n < would be problematic for integrals on < r < oo involving e~ nr dr, 
but these integrands are replaced by ones involving e~ xr dr with x > 1. 

5. Asymptotics OF arg J„ + iY n AND of the zeros of Bessel functions 

Denote H n (x) — J n {x) + iY n {x). It is well known that for all real numbers n, 
\H n (0 + )\ — oo and \H n (x)\ is a decreasing function of x for x > (Watson [W] p. 
446). Moreover, 

i^J-lm^ - Jn{x)2 + Yn{x)2 - nx]Hn{x){2 > 

Therefore, for all real values of n, as x > increases, H n (x) traces a simple spiral 
counterclockwise in the complex plane. To choose a well-defined branch of the 
argument 

9(x,n) = a,ig(H n (x)) 

note that Y n (x) = ( J n (x) cos(nn) — J_„(a;))/ sin(ri7r) and one has the asymptotic 
formula J n (x) — (x/2) n /T(n + 1) as x -> 0+. We deduce that for n > 0, 

H n (x)/\H n (x)\ ->•-* = er ln/2 as x 0+ 
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whereas for < n < —1/2, 

H n (x)/\H n {x)\ -> - sin nn-i cos nn = e -Kn+l/2)ir as x Q + 

Therefore, for x > sufficiently small H n (x) is in the 4th quadrant, and a consistent 
definition of the branch is given by 

0(0+, ra) = -tt/2 (n>0) 

and 

6>(0+,n) = -(n+ 1/2)tt > n > -1/2 
It then follows that the mth positive zero of J n (x) satishes 

(5.1) 9(p(m,n),n) — run — n/2 

for m = 1,2,.... This implicit equation can then be used to extend the definition 
of p{m, n) for all real m and n satisfying m + n > 0. 

Next, we establish a few preliminary upper and lower bounds for p(m,n). It 
follows from the fact that d x 9(x,n) > that p(m,n) well-defined and infinitely 
differentiable. In the range |n| < 1/4, the formula for d x 9(x,n) and the estimates 
\H n (x)\ ■< a; -1 / 4 in < x < 1 and \H n (x)\ < x~ x l 2 in 1 < x < oo imply that 

(5.2) p(m,n)^m (m > 1/2, |n| < 1/4) 

In the range n > 1/4, m > 3/4, we prove (well-known) upper and lower bounds 

(5.3) p(m,n) — n ~ m + m 2 ! 3 ^! 3 

Let n > 0, then according to [W] (p. 485-487), n < y n where y n = p(l/2,n) is the 
smallest positive zero of Y n (x). Since 9{y n , n) = 0, 

np{m,n) 

mir — 7r/2 = / d x 9(x,n)dx 



For all n > 1/4, |if„(a;)| r^ 1 / 3 in n. < a; < n+ n 1 / 3 and |iT n (x)| x^ 1/4 (x - 
n)" 1 / 4 in i > n + n 1 / 3 . These estimates along with the formula for d x 9(x,n) above 
yield (|S3]) . 

To prove Theorem l2.31 we require symbol estimates for 0(x, n) expressed in terms 
of the variables {f3,v) with cc = i^sec/3, n = v. 

Lemma 5.1. Denote 

a(v, ff) = 9(usec0, v) 

There is an absolute constant C such that if Cv~ x / 3 < j3 < tt/4, then 
a) \did k p [a{P, v) - z/(tan/3 - p)]\ ^ i/-J/3 -fc j + fc > 1 

\did k {a(p,v)\ <v l - j ^- k j + k>l 

In the transition region, where c < vf3 3 < C, for some absolute constants < 
c < C < oo more detailed asymptotics are required. Denote y = [y 2 ! 3 tan 2 /3)/2 
and define the derivative of arg A(y) by 

v ' My) \A{y)\ 2 

where A(y) — u(y) + iv(y) is the Airy function in (|4.13[) . The asymptotic expansion 
for a(v, j3) is given by 
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Lemma 5.2. If j + k > 1, then 

a) d v u = l(v~ 1/3 tan 2 P)B{y){\ + 0{p 2 + u~ 4 / 3 p~ 2 )) 

b) d a = {v 2/3 tirn /3 sec 2 f3)B{y){l+0(f3 2 + v- 1 /3- 1 )) 

Proof of Lemma 15.11 First note that the definition of a and (|4.1I) imply 

c»,3a= (sec 2 p- l)u + Im^ 

and 

d v o = (tan p-fl+hn — 

a 

Thus Lemma 15.11 follows from Proposition 14.11 and Corollary 14.31 
Next, we prove Lemma 15.21 

d v G = tan — + lm(d v a/a) 

= tan ,3 - ft + lm.(d v b/b) + lm[{d v a - d v b)/a] + lm[(d u b)(l/a - 1/6)] 
= tan/3 - + lm(8 u b/b) + 0{\d v (a - b)\/\a\ + \d n b\\a - b\/\ab\) 
= t&np-p + \m{d v b/b) + 0((p + z/" 1 / 3 ) 5 ) 
Moreover (with y = {v 2 ! 3 tan 2 P)/2), and recalling that B(y) = Jm(A'(y)/A(y)), 

Im = - tan 3 /3/3 + ^— ^- tan 2 pB(y) 
b 3 

and tan/3 - P - (tan 3 /3)/3 < P 5 . Thus, 

= tan 2 /3B(y) + 0((/3 + ^ 1 / 3 ) 5 ) 

Similarly, 

c^er = t/(sec 2 /3 — 1) + lm.{dpa/a) 

= iv[scc 2 P-1 + lm(dpb/b) + Im[(<V - dpb)/a] + kn[(d b)(l/a - 1 /&)] 
= ^[sec 2 /3 - 1] + Im(^6/6) + 0(i/(/3 + i^ 1/3 ) 4 ) 
Moreover, 

Im-2- = -^tan 2 /3sec 2 ^ + i/ 2/3 (tan/3)(sec 2 /3)B(y) 
6 

Since i/[sec 2 /3 — 1] — z/tan 2 /3sec 2 /3 ^ i^/3 4 , 

<V = i; 2 / 3 (tan /3) (sec 2 /8)fl(y) + 0(y{p + z^ 1/3 ) 4 ) 

Finally, in order to write the error as a multiplicative expression, we use |-B(y)| ~ 
(1 + y) 1 / 2 which follows from the upper and lower bounds on A{y) and A'(y) (and 
ultimately from the Wronskian formula). 

We can now deduce Theorem 12.31 from Lemmas 15.11 and 15.21 Consider the func- 
tions v and p of n and to defined by 

v = n; cos P = nj p{m, n) 

For 3/4 < m < 3n, p(m,n) — n w to 2 / 3 ^ 1 / 3 implies 

P(m,n) w (m/n) 1/3 if 3/4 < m < 3n 

i//3 2 « m 2/ V /3 /3w(m/n) 1/3 
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In particular when m > 3/4, (3 > cn -1 / 3 for some small absolute constant c > 0. 
We will be treating three ranges of (3 in different ways. One is the transition region 
where cnT 1 ! 3 < (3 < Cn -1 / 3 where C is a large constant. The second is the region 
Cn- 1 / 3 < p < p for some small absolute constant j3o, and the third is the classical 
region (3o < (3. We will never need to consider (3 smaller than cn~ 1//3 . 
Differentiating a{v, f3) — mit — 7r/2, we find the implicit formulas 

o n p = -5 ; o m /3 = - — 

dpa dpa 

We are going to prove by induction the property Q(j, k) that says that for all 
< £1 < j and all < £ 2 < k 7 

(/? « (to/ti) 1 / 3 ) The property Q(0, 0) is trivial. For Q(0, 1), observe that by Lemma 



d/3 a Zv z l 6 tan p sec (3 

Moreover, 

The proof of Q(1,0) is similar. Suppose that Q(j,k) is valid. Differentiating the 
implicit formula for 8-^8^(3 with respect to m, there are three types of things that 
can happen. First the derivative falls on the denominator, which is a power of dpa, 
in which case the expression is multiplied by a constant times 

(d%a)d m p 
dps 

But recall that \dpa\ >z P 2 v, \d^a\ ^ (3v, and \d m f3\ ^ f3~ 2 i>~ 1 so that the product 
is majorized by f3^ 3 i/~ 1 . This is the new factor required for the estimate Q(J + 1, k). 
If the derivative falls on the numerator, then it may increase the degree of differen- 
tiation on a derivative of /3, but always below the level of the induction hypothesis. 
Replacing a derivative of /? by one derivative higher yields a change in estimation 
of the whole by the appropriate factor f3~ 3 v~ l . Finally, the differentiation may 
land on a derivative of Da. This replaces Da by (dpDa)d m (3, so the change in 
the estimation is the same (difference between Da and dpDa is a factor f3^ 1 and 
d m f3 is bounded by ft^v^ 1 . Again the product is fi~ z v~ x , which is the factor we 
want. Similarly, to prove Q(j, k + 1), differentiation with respect to n produces an 
estimate that differs by the factor v~ x from the bound for Q(j, k). 

Next we can use Q(j, k) to prove the property P(j, k) that for all < £\ < j and 
all < £ 2 < k, 

\d^(p(m,n) - n)\ £ /3 2 - 3 V"^ 2 (~ m^^n 1 ^) 

In fact, implicit differentiation of pcos/3 = n with respect to n and with respect to 
m gives 

d n p — 1 = (sec f3 — 1) + p(tan (3)d n (3 

and 

d m p = p{t'&n (3)d m (3 

The induction argument is similar to the proof of Q(j, k) and is left to the reader. 
This concludes Theorem 12.31 (a). 

21 



We turn now to the lower bound (Theorem 12.31 b) . 



d n p — 1 = (sec f3 — 1) + p(tan f3)d n /3 

d o 

= (sec j3 — 1) — p(tan/3) — — 

= + 0(/? 4 ) - i sin (3 tan 0(1 + 0(/3 2 + i/~ 4 / 3 /T 2 ) 

1^+0(^4^-4/3) 

6 

Note that we need the precise asymptotics because the lower bound of order 
j3 2 fa m 2 / 3 n~ 2 / 3 requires the coefficient 1/2-1/3 = 1/6 > 0. The error term 
is lower order if j3 « 1 and n~ 4 / 3 << /3 2 . Thus we have Theorem 12.31 (b) for n 
sufficiently large, in which case i/~ 2 / 3 << j3 < CV -1 / 3 . We do not need smaller 
values of (3, because for m > 1/20, and j3 defined implicitly in terms of p(m,ri), 
/?« (m/n) 1 / 3 » is- 2 ' 3 . 

At last, here are some details of the much simpler estimates in what we are 
calling the classical region. Let n > 1 and x > (1 + c)n. Then we are in the range 
in which (x — n) ~ x. As suggested in the remark of [D] p. 1069, following the 
methods above one finds 

H n (x) ^ x- 1/4 (x - n)- 1/4 e lixsinp - n0) s{x,n) 

where fi(x, n) is defined by xcos/3 = n and s(x,n) satisfies 

\did*s(x,n)\<x-i- k 

Moreover, it is well-known that |s(x,n)| > c > 0. It follows that 6 (x, n) — axgH n {x) 
(defined using any appropriate branch) satisfies 

9(x, n) = x sin j5 — n(3 + E(x,n) 

where E(x, n) is a symbol satisfying 

\did k n E{x,n)\<x-i- k 

It is not hard to extend this estimate to negative values of n. Indeed the remaining 
range of n is the range —1/2 < n < 1. For a fixed range of the parameter n, 
the Hankel formula for H n (x) (jW] pp. 196-198) may be used. Let x > 1 and 
— 1/2 < n < xcospo, where /3q is a small, fixed constant. Then since po > 0, 
x > (1 + c)n for some c > and we are in the range in which (x — n) ~ x. 
One calculates that 

= 9 U a ; d n j3 



xsin/3 
and 

d n e = -(3 + d n E- : d x e + d x E 

Differentiating the implicit equation for p(m, n) 

9{p(m 1 n), n) — mir — n/2 
with respect to m and n, we find 

7T 7T 
dmP = 5 ^7 p x T = -r—B + F l(P, n ) 

o x t){p[m,n),n) smp 
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and 

g _ d n 9(p(m 7 n),n) _ (3 + 

d x 8(p(m,n),n) sin/3 ' 

with F\ and F% satisfying 

and pcos/3 = n, /3 </3<7r/2 + e(/3 extends a small amount beyond 7r/2 to 
accomodate the negative values of n — we only care about p sufficiently large, so 
e can be arbitrarily small.) 

Finally, this implicit asymptotic formula for p(m, n) can be differentiated many 
times expressing derivatives of pirn, n) in terms of lower derivatives. This yields 

\d> m d*p(m,n)\3p 1 -S- k Km L -S- k 

for all m > 100 + cn for a small fixed c > 0. For the asymptotic gradient we already 
have the implicit formula 

(d m p, d n p) = (vr/ sin/3, (3/ sin/3) + F(p, n) 

with F a symbol of order -1 as above. Thus of order 0(1 /m) since p ss m in this 
range. We also have the implicit equation for p that gives us /3 as a function of p 
as follows: 

n(tan/3 - /3) + E(p, n) = mn - tt/2 

so that (p cos /3 = n) 

Tfl IT Tfl IT 

tan/3 - p = (1-(E + 7t/2)/7tto) = (1 - 0(1/ m)) 

n n 

Thus if we define a(m, n) by 

tana — a — irm/n 

(across tt/2 as needed is perfectly ok) then a = f3 + 0(1 /m) and we get the as- 
ymptotic formula we wanted. (We can also characterize the error term as a symbol 
rather than with a bound, but we don't need this.) 

6. Final Remarks 

The length spectrum of the disk is far from generic. The regular fc-gon trajectory 
on the disk can be rotated around the circle giving a one-parameter family of 
periodic geodesies of the same length, whereas in general the length of the geodesic 
with k reflections varies depending on where its vertices are. Nevertheless, we expect 
that the contribution to the wave trace at times t greater or equal to the perimeter 
from this family of fc-reflection periodic trajectories resembles the effect of the family 
of regular fc-gons in the case of the circle. Indeed the almost integrable behavior of 
the dyamical system of geodesic flow on convex planar regions near the generalized 
geodesic that follows the boundary has been treated in detail by Lazutkin [L] and by 
Melrose and Marvisi [MM]. Moreover, as we saw above, the asymptotics depended 
in a fundamental way on nondegeneracy conditions on derivatives of Airy functions. 
This is a hopeful sign, since general microlocal constructions of the parametrix of 
the wave equation near the boundary also involve Airy functions |AM[ IMS! IMTj . 

As already mentioned, certain symbol properties of Bessel functions J v (x) as a 
function of the two variables x and v were already proved in Ionescu and Jerison in 
[L7] . There are three features of the treatment in 1J that are not enough for our 
purposes here. First, there is the minor point that we need bounds on derivatives of 
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all orders, not just the first two. Second, the full symbol-type estimates (conjectured 
in [D], P- 1069) while valid are not sufficient here. Those bounds are for derivatives 
in the variables (x,v), whereas for the present purpose in the range v < x < 2v, 
especially when x is near v, it is necessary to distinguish a special directional 
derivative, namely, the derivative with respect to v with the ratio x/v held fixed 
(or equivalently with j3 defined by cos/3 = v/x held fixed). Third, [LJ] treats only 
upper bounds on J v {x) +iY v (x) and its derivatives. We need both upper and lower 
bounds on the argument and its derivatives, which require detailed asymptotics, 
not just bounds on the Bessel functions and their first derivatives. 

References 

[AM] K. G. Andersson and R. B. Melrose, The propagation of singularities along gliding rays, 
Invent. Math. 41 (1977) no. 3, 97-232. 

[C] J. Chazarain, Construction de la parametrix du probleme mixte hyperbolique pour 
l'equation des ondcs, C. R. Acad. Set. Paris Ser. A-B 276 (1973), A1213-1215. 

[CdV] Y. Colin de Verdiere, Spectre du laplacien et longueurs des geodesiques periodiqucs II, 
Compositio Math. 27 (1973), 159-184. 

[DC] J. J. Duistermaat and V. Guillemin, The spectrum of positive elliptic operators and periodic 
bicharacteristics, Invent. Math. 29 (1975) 39-79. 

[GM] V. Guillemin and R. B. Melrose, The Poisson Summation formula for manifolds with bound- 
ary, Adv. in Math. 32 (1979) 204-232. 

[HeZ] H. Hezari and S. Zelditch, Inverse spectralproblem for analytic (X/2Z) n -symmetric domains 
in R™, arXiv: 0902.1373 

[IJ] A. D. Ionescu and D. Jerison, On the absence of positive eigenvalues of Schrodinger oper- 
ators with rough potentials, Geom. Fund. Anal. 13 (2003), no. 5, 10291081. 

[L] V. F. Lazutkin, Construction of an asymptotic series of eigenfunctions of "bouncing ball" 
type, Proc. Steklov Inst. Math. (1968) 125-140. 

[MM] S. Marvisi and R. B. Melrose, Spectral invariants of convex planar regions, J. Diff. Geom. 
17 (1982) 475-502. 

[MS] R. B. Melrose and J. Sjostrand, Singularties of boundary value problems, Comm. Pure 
Appl. Math. 31 (1978) 593-617. 

[MT] R. B. Melrose and M. Taylor, Near peak scattering and the corrected Kirchkoff approxi- 
mation for a convex obstacle, Adv. in Math. 55 (1985) 242-315. 

[W] G. N. Watson, A treatise on the theory of Bessel functions, 2nd edition, Cambrige U. Press 
(1966) 

[Z] S. Zelditch, Inverse spectral problems for analytic domains II: domains with one symmetry, 
Ann. of Math. (2) 170 (2009) 205-269. 

Institut de Fourier. U. de Grenoble I 

E-mail address: yves.colin-de-verdiereOujf-grenoble.fr 

Massachusetts Institute of Technology 
E-mail address: vwg0math.mit.edu 

Massachusetts Institute of Technology 
E-mail address: jerisonOmath.mit.edu 



24 



